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ABSTRACT 

We present spectroscopic observations of red giant branch (RGB) stars over a large expanse in the 
halo of the Andromeda spiral galaxy (M31), acquired with the DEIMOS instrument on the Keck II 
10-m telescope. Using a combination of five photometric/spectroscopic diagnostics — (1) radial veloc- 
ity, (2) intermediate-width DD051 photometry, (3) Nai equivalent width (surface gravity sensitive), 
(4) position in the color-magnitude diagram, and (5) comparison between photometric and spectro- 
scopic [Fe/H] estimates — we isolate over 250 bona fide M31 bulge and halo RGB stars located in 
twelve fields ranging from R — 12-165 kpc from the center of M31 (47 of these stars are halo members 
with R > 60 kpc). We derive the metallicity distribution function of M31 RGB stars in each of these 
fields by comparing the stellar location in the (I, V — I) color-magnitude diagram to a finely spaced 
grid of theoretical isochrones. The mean of the resulting M31 spheroid (bulge and halo) metallic- 
ity distribution is found to be systematically more metal-poor with increasing radius, shifting from 
([Fe/H]} = -0.47±0.03 (a = 0.39) at R < 20 kpc to ([Fe/H]) = -0.94±0.06 (a = 0.60) at R ~ 30 kpc 
to ([Fe/H]) = -1.26±0.10 {a = 0.72) at R > 60 kpc, assuming [a/Fe] = 0.0. These results indicate 
the presence of a metal-poor RGB population at large radial distances out to at least R = 160 kpc, 
thereby supporting our recent discovery of a stellar halo in M31: its halo and bulge (defined as the 
structural components with R~ 2 power law and de Vaucouleurs i? 1 / 4 law surface brightness profiles, 
respectively) are shown to have distinct metallicity distributions. If we assume an a-enhancement 
of [a/Fe] = +0.3 for M31's halo, we derive ([Fe/H]) = -1.5±0.1 [a = 0.7). Therefore, the mean 
metallicity and metallicity spread of this newly found remote M31 RGB population are similar to 
those of the Milky Way halo. 

Subject headings: galaxies: individual (M31) - galaxies: structure - Galaxy: abundances - techniques: 
spectroscopic 



1. INTRODUCTION 

Large galaxies such as the Milky Way and Andromeda 
(M31) are believed to ha ve been assembled hierarchi- 
cally ijSearle & ZinrJlT978() . The growth of such galax- 
ies is powered by the continual accretion of smaller 
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dwarf galaxies that a re tidally destroyed as the y fall into 
the larger potential ijZentner & Bullockll2003tL Numer- 
ical simulations suggest that the most massive merging 
events occur early on (^lookback ^ 8 Gyrs) and form the 
inner halo (R < 20 kpc), whereas t he recently accreted 
dwar fs form structure in the halo ijBullock fc JohnstorJ 
2005). This formation scenario leads to several predic- 
tions that can be observationally tested. For example, 
accretion of dwarf galaxies should naturally produce an 
extended stellar halo in massive hosts. Furthermore, the 
recent infall and subsequent tidal disruption of dwarf 
satellites should produce a large amount of stellar sub- 
structure (i.e., tidal streams) still existing within galaxy 
halos. Finally, since the most massive merging events 
that formed the inner parts of hosts li ke the Milky Way 
and M31 are also the most metal-ric h llFont et alJl200fit 
[ Robertson et alJ 120051 iR.enda et all 120051 iBrook et alJ 
2004), this formation scenario suggests that the inner 
parts of massive galaxies should be chemically different 
from their halos. 

Recently, a large number of ground and space based 
observations have targeted both the Milky Way and 
M31 to try to confirm these predictions directly. In 
this respect, M31 is often a better choice than our own 
Galaxy since we have a global, external view of it. The 
wealth of data collected to date has indeed confirmed cer- 
tain theoretical predictions, while, at the same time, it 
has raised some new puzzles. For example, wide-field 
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Fig. 1. — Location of the fields used in this study relative to M31 (center). The size of M31's visible disk, and orientation, is illustrated with 
two solid lines in the center. Lines of constant longitude and latitude arc also shown. The smaller gridded squares represent schematically 
the field centers of the KPNO fields for which we have DD051 photometry and the two larger squares indicate the approximate positions 
of the CFHT fields without DD051 photometry. The filled circles show the approximate locations of the twelve fields for which we have 
Keck/DEIMOS spectroscopic observations. 



star-count maps of both galaxies have revealed abun- 
dant substructure in the halos of these systems, as pre- 
dicted by s imulations. In the Milky Way, the Sag ittar- 
ius stream (llbata et alJl!994t iMaiewski et aUl2003ft. the 
Magellanic stream IjMathews on^^lgarVjfc M urravll974l) 
and the Monoceros strearn f^mnvet^nTl2003() present 
clear evidence that merging processes are prominent in 
massive galaxies and are still oc curing today. Simi- 
larly, star-count m aps ofM31 fe.g-. lFerguson et aDl2002t 
llbata et al.ll2005|) have shown that its spheroid is very 
imhomogenous with many prominent densi ty enhance- 
ment s (e.g., the giant southern stream — llbata et al.l 
l200l . 

The comparison of the Milky Way and M31 halo has 
also revealed several unanticipated results and notable 
differences between the two galaxies. Despite their simi- 
lar overall size, the stellar density of M31's spheroid ap- 
pears to be 10 x higher than that at a comparable loca- 
tion in the Milky Way fReitzfiL Guhatha kurta. & Gouldl 



11998ft . iMould fc Kristianl {1986) surprisingly found that 
the "halo" of M31 at R ~ 7 kpc is also 10 x more 
metal rich t han the Milky Way's. This w as also 
measured bv iDu rrell. H arris, fc Pritchetl l)1994(l and by 
iRich . Mighe ll. fc Neill <|1996l) over a larger radial dis- 
tance, from R = 5.3-19.4 kpc. Further studies of the 
metallicity distribution function (MDF) in M31 not only 
confirmed this result, but also showed that there is no ev- 
idence for an abundance gradient out to 30 kpc (Durrell, 
Harris, & Pritchet 2001, 2004, Bellazzini et al. 2003), 
a result at odds with model predictions. Differences be- 
tween the Milky Way and M31 are also seen in the distri- 
bution of the ages of stars at large radii. The canonical 
picture of an old, metal-poor popul ation simply d o es not 
appear to describe M31's "halo". iBrown et aD (|2003f) 
used the Advanced Camera for Surveys (ACS) on the 
Hubble Space Telescope (HST) to target a minor-axis 
field at R ~ 10 kpc in M31 and produced a superb 
color-magnitude diagram (CMD) extending well below 
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the main-sequence turnoff. Their data suggest a broad 
distribution of stellar ages in M31's "halo": over half 
the population of stars in this field have ages < 11 Gyr 
old. By mass, 30% of the stars are found to be 6-8 Gyr 
old. Rich et al. (2006, in preparation) have now ob- 
tained spectroscopy of stars in this HST field using the 
DEIMOS spectrograph on the Keck II telescope. We use 
these data below to confirm that the RGB population in 
the Brown et al. M31 "halo" field is much more metal 
rich than RGB stars in the Milky Way halo. 

Until recently, there appeared to be another striking 
structural difference between the Milky Way and M31 
spheroids. W hereas the Milky Way shows a de Vau- 
couleurs i? 1 / 4 (Ide Vaucouleursl H958) surface brightness 
profile in the inner regions of the Galaxy and a power 
law R~ 2 projected profile for the halo in the outer parts, 
iPritchet fc van den Berghl l)1994l) found that the entire 
spheroid of M31 could be fit by a single i? 1 / 4 profile. 
Their star-count measurements extended out to R ~ 20 
kpc along the south-east minor axis of M31, and demon- 
strated that the surface brighness profile of M31 falls 
off very steeply with increa sing radius near the limit 
of the ir survey. Similarly, iDurrell. Harris, fc PritcheH 
( 2004) found that an i? 1 / 4 law is consistent with M31's 
s urface brightness profile out to 30 kpc. 

Ost heimerl l)2003|) carried out a survey that was a signif- 
icant improvement over previous work in terms of both 
spatial coverage (R ~ 10-165 kpc), photometric depth 
(1.5 mag arcsec -2 fainter), and signal-to-background. 
His M31 surface brightness profile showed the first sign of 
a flattening (e.g., deviation of the slope from a de Vau- 
couleurs profile) in t he outermo s t few bins (>80 kpc). 
By following up the lOstheinierl (2003) imaging obser- 
vations with multiobject Keck/DEIMOS spectroscopy, 
we have now been able to verify M31 red giant branch 
(RGB) member stars in each of the Ostheimer fields out 
to R ~ 160 kpc. The spectroscopic data unequivocally 
show a break in the surface brightness profile of M31. In 
iGuhathakurta et alJ l|2006aO . we present the first detec- 
tion of an R~ 2 halo of stars in M31 extending out to a 
projected radius >160 kp c. Based solely on photometric 
data out to R ~ 55 kpc, llrwin et all l|2005D also find a 
break in the surface brightness profile of M31 (although 
they claim that the outer component does not resemble 
a population II halo). 

Following up on this discovery, this paper presents the 
first detection of a metallicity gradient in M31 9 . Further- 
more, we show that the crossover between the metal- 
rich and metal-poor components in our sample occurs 
at a minor-axis distance of ~30 kpc, in excellent agree- 
ment with the transition radius that separates the newly 
discovered halo of M31 from the inner de Vaucouleurs 
(i? 1 / 4 ) bulge 10 . We propose that this new component is 
in fact the stellar halo of M31 (see also Chapman et al. 
2006). Taken together, these results clarify that the ac- 
tual structural disparities between the Milky Way and 

9 Although Ostheimer 1 2003) did not claim to have seen an obvi- 
ous/significant metallicity gradient in M31, he did provide evidence 
for an increase in the fraction of metal-poor stars in his outermost 
annuli. 

10 Throughout the rest of this paper, we use the terms "bulge" 
and "halo" to refer to M31's i? 1 / 4 and R~ 2 structural components 
which dominate the R < 30 kpc and R > 30 kpc regions of the 
galaxy, respectively. 
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Fig. 2. — Sample M31 RGB spectra, one star from each of 
our twelve fields, as well as the spectrum of a typical foreground 
Milky Way dwarf star (bottom). The spectra have been corrected 
to zero velocity, normalized at A ~ 8500A, offset in y arbitrarily, 
and smoothed using a 10 pixel (3A) boxcar function for illustration 
purposes only. The locations of the Call triplet: 8498, 8542, and 
8662A and pressure sensitive Nal doublet: 8183 and 8195A are 
indicated with arrows at the top of the figure. As expected, the 
Nal doublet is measured to be very strong in cool foreground Milky 
Way dwarf stars, but is absent in M31 RGB stars. 



M31 may not lie in the properties of their halo popu- 
lations, but rather in the relative sizes of the bulges of 
the two systems. The discovery of this metal-poor stel- 
lar halo also provides a powerful confirmation of galaxy 
formation models (see § I7.5[) . 

2. OBSERVATIONS AND DATA REDUCTION 

The locations of the imaging/spectroscopic fields pre- 
sented in this study are shown in Figure 1. The wide- 
field imaging observations were obtained with the Kitt 
Peak National Observatory (KPNO) 4-m telescope an d 
Mosaic camera for most of our fields (|Ostheimerll2"003]l . 
These are shown as gridded squares. The two larger open 
squares show the positions of very wide- field imaging ob- 
servations obtained with the Canada-France-Hawaii 3.6- 
meter telescope and MegaCam camera. The filled circles 
indicate the positions of our twelve spectroscopic fields, 
all obtained with the Keck II telescope and DEIMOS 
instrumen t. The imaging fo r the outer fields was ob- 
tained by lOstheimerl l|2003(l using the KPNO Mosaic 
camera, in the Washington System M a nd Ti bands and 
the i ntermediate-width DD051 filter (jMaiewski et alJ 
2000). This filter allows measurement of the surface- 
gravity sensitive Mgb and MgH stellar absorption fea- 
tures and therefore provides a means to discriminate fore- 
ground Milky Way contaminants from M31 RGB stars. 
The innermost two fields in our study, Hll and HI 3s, 
are near enough to M31 that there is minimal foreground 
Milky Way dwarf contamination. We therefore did not 
observe these fields in the DDObl filter. Rather, we ob- 
tained very wide-field Canada-France-Hawaii Telescope 
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TABLE 1 



Date 


Mask 


Pointririj 


y center! 


Field PA 


No. Sci. 






a J2000 


^J2000 


(°E of N) Targets 1 






(h .m .s \ 
V ■ • J 


IO.I.IIS 

\ ■ ■ J 






2002 Ane- ItS 


a3_l 




-1-39-02-39 2 


+64.2 


85 


2(102 Alio- Ifi 


a0_l 


00-51-51 32 


4-39-50-21 4 


-17.9 


89 


2(102 Ort 1 1 


a3 2 


00-47-47 94 


Toy . uo . tju . o 


+178.2 


80 


2002 Ort 12 


a0_2 


00-51 -29 59 


T J3 ■ T:^t. O 


+90.0 


89 


2003 Son 30 


al3 1 


00-49- ^4 


tou . Oij . iy . o 


+0.0 


80 


2003 Spn 30 


al3_2 


00-41 -28 27 




+0.0 


71 


2003 Son 30 


mll_l 


01 -29-34 44 


4-34-1 3-45 4 


+0.0 


72 


2003 Ort 1 


mll_2 


01 -29-^4 35 


4-34-27-45 5 


+0.0 


68 


2003 Oct 1 




01-09-51 75 


4-37-46-59 8 


+0.0 


75 


2003 Ort 2fi 


a3_3 


00-48-23 1 7 


4-39-1 2-38 5 


+270.0 


83 


2004 Tune 1 7 

ijVJU^t U LlllL _1_ I 


a0_3 


00-51 -50 4fi 


4-40 07 00 9 


+0.0 


90 


2004 Spn 20 


Hll_l 


00-46-21 02 


4-40-41 -31 3 


+21.0 


139 


2004 Sep 20 


Hll_2 


00:46:21.02 


4-40-41 -31 3 

I t: \J i t: 1 i 'J 1 . 'J 


—21.0 


138 


2004 Sep 20 


H13s_l 


00-44-14 7fi 

UW-T:T:t It:- 1 U 


4-39-44-18 2 


+21.0 


134 


2004 Sep 20 


H13s_2 


00-44-14 7fi 


4-39-44-18 2 

| ^ J .'1'l.lUtiJ 


—21.0 


138 


2005 Jun 9 


m6_2 


01:08:36.22 


4-37-28-59 6 


+0.0 


72 


2005 Jul 7 


m8_l 


01:18:11.56 


+36:16:24.9 


+0.0 


56 


2005 Jul 7 


m8_2 


01:18:35.87 


+36:14:30.9 


+0.0 


59 


2005 Jul 8 


mll_3 


01:30:01.53 


+34:13:45.4 


+0.0 


80 


2005 Jul 8 


mll_4 


01:30:37.33 


+34:13:27.4 


+0.0 


75 


2005 Aug 29 


al9.1 


00:38:16.05 


+35:28:07.2 


-90.0 


71 


2005 Sep 6 


d2_l 


01:17:07.46 


+33:29:25.1 


-90.0 


139 


2005 Sep 6 


d2_2 


01:16:43.29 


+33:34:25.8 


+0.0 


111 


2005 Sep 7 


bl5.1 


00:53:23.63 


+34:37:16.0 


-90.0 


65 


2005 Sep 7 


bl5.3 


00:53:37.77 


+34:50:04.1 


-90.0 


74 


2005 Sep 8 


d3.1 


00:36:03.83 


+36:27:27.4 


+90.0 


120 


2005 Sep 8 


d3.2 


00:35:39.61 


+36:21:41.8 


+0.0 


122 



1 Some targets were observed on multiple masks. 



(CFHT) photometry in the <?' and i' filters using the 
1 degree 2 MegaCam camera. The g 1 and i' photometry 
was converted to Johnson-Cousins V and / using stan- 
dard Landolt field observations. 

Details on the imaging observations, slitmask de- 
sign, spectroscopic observations, and data reduction of 
the fields aO, a3, al3, d3, a!9, m6, b! 5 , m8 , d2, 
and mil are given in § 2 of iGilbert et alJ 1)20061) and 
iGuhathakurta et alJ| 2006b). Similar information for the 
two innermos t fields Hll and H13s can be found in § 3 of 
iKalirai et all <|2006f) . Keck/DEIMOS spectroscopic ob- 
servations were obtained in each of t he ab ove pointings 
as dis cussed in IGilbert et alJ (|2006) and IKalirai et al.1 
(2006). We briefly summarize the information presented 
in these papers. We used the 1200 lines mm -1 grating 
(dispersion = 0.33A pixel -1 ), providing a spectral resolu- 
tion of 1.3A (FWHM) for typical 0'/8 FWHM seeing. We 
targeted the brightest M31 RGB stars (20.0 < I < 22.5 
- the brightest in this window may in fact be asymptotic 
giant branch stars) in the photometry and built masks 
in an iterative process that maximized the total num- 
ber of slits selected for highest priority objects. These, 
in the case of the outer fields, were largely based on 
the position of the star in the CMD and the position 
in the (M - DD051) versus (M - T 2 ) color-color dia- 
gram (e.g., Majewski et al. 2000; Palma et al. 2003). 
As mentioned earlier, the latter provides a preselection 
to ensure high probability RGB stars, which is especially 



important in the outer regions of M31 where true RGB 
stars have sparse density. For the inner fields, a combi- 
nation of the magnitudes and a stellarity (star-like) limit 
based on the morphology of the sources from SExtractor 
l)Bertin fc Arnoutslll996ri was used to prepare the masks. 
In general, the high density of objects in the inner fields 
allowed us to target ~150 targets on each mask. All of 
the spectra were inspected visually and assigned a quality 
control index based on the quality and number of abso rp- 
tion lines visible. As discussed in IKalirai et al.l (J2006) , a 
significant number of the fainter targets yielded spectra 
that did not show any obvious features due to low signal- 
to-noise (S/N). We determine the velocities for all objects 
that contain at least two spectral features (at least one 
definite and one other marginal line) by cross-correlating 
the observed spectra with a large database of template 
stellar and emission- and absorption-line galaxy spectra. 
The mean S/N of these good spectra is ^10 per pixel 
and the velocity uncertainty from the cross-correlation 
is empirically estimated to be ~15 km s . The mean 
magnitude for these stars (with S/N ~ 10 per pixel) is 
I - 21.2. 

After removing galaxies, the twelve fields (27 masks) 
in this study contain 1070 stars for which we obtained a 
reliable velocity. Figure [2] shows a single representative 
spectrum of an M31 RGB star from each of these fields. 
Table 1 presents a summary of the observations. 

3. A CLEAN SAMPLE OF M31 RGB STARS 

A fraction of the 1070 stars in our sample are in fact 
foreground Milky Way dwarfs. In order to measure the 
MDF in M31's bulge and halo, we first need to isolate 
M31 RGB stars from this contamination. We have de- 
veloped a sensitive technique that uses probability distri- 
bution functions (PDFs) calculated from a training set 
of kno wn RGB and dwarf stars to provide this discrim- 
inant ijGilbert et alJl2006() . Using five criteria — (1) ra- 
dial velocity, (2) intermediate-width DD051 photome- 
try, (3) Na I equivalent width (surface gravity sensitive) , 
(4) position in the CMD, and (5) comparison between 
photometric and spectroscopic [Fe/H] estimates — each 
star in our sample is assigned five likelihood values of be- 
ing a giant or a dwarf based on its location within each 
diagnostic plot. The individual probabilities are then 
combined to yield the final discriminant of whether the 
star is truly an M31 RGB star. In Figure |3| we show the 
properties of confirmed M31 RGB stars in six of our outer 
spectroscopic fields, al3, al9, m6, bl5, m8, and mil. We 
have plotted four of the five diagnostics used to distin- 
guish RGB stars from dwarfs, and have also overlayed 
empirical PDFs derived from the training set of definite 
RGB (solid darker curve) and dwarf (dashed/thin curve) 
stars. It is clear that the outer RGB stars predominantly 
agree with the RGB training set in all of their properties. 

A detailed discussion of the procedure used to average 
th e individual PDFs in to a likelihood distribution is given 
in IGilbert et alJ ((2006) . In Figure 0] we present the final 
distribution of weighted average likelihoods, (Li). Al- 
though any star with (Li) > is a preferred M31 RGB 
star (whereas stars with (Li) < are likely to be Milky 
Way dwarfs), we only select stars with (Li) > 0.5 for 
our clean M31 RGB sample. This ensures that the star 
is three times more likely to be an M31 RGB member 
than a Milky Way dwarf (see Gilbert et al. 2006). We 
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Fig. 3. — Properties of confirmed M31 RGB stars in six of our outer spectroscopic fields (al3, al9, m6, bl5, m8, and mil - all located at 
R > 60 kpc). Four of the five diagnostics used to distinguish RGB stars from foreground Galactic dwarf star contaminants are illustrated, 
(a) Radial velocity distribution of M31 RGB stars and dwarf stars [shaded-bold and dot-dashed histograms, respectively], compared to 
empirical probability distribution functions (PDFs) derived from training sets consisting of definite RGB and dwarf stars in all of our fields 
(solid and dashed curves, respectively), (b) Same as (a) for the DD051 parameter, which is based on the star's location in (M — DD051) 
vs. (M — T2) color-color space. For both (a) and (b), we have scaled back the dwarf star histograms by a factor of three for clarity and 
arbitrarily scaled the training set PDFs. (c) Equivalent width of the Nal 8190A absorption band versus de-reddencd (V — I)o color for the 
RGB and dwarf stars [bold circles and crosses, respectively]. Bold and thin contours show 10%, 50%, and 90% e nclos ed fractions for RGB 
and dwarf star PDFs, respectively, (d) Same as (c) for the photometric metallicity estimate (CMD based, see S I5.1I versus spectroscopic 
metallicity estimate (based on the strength of the 8500A Call triplet, see § 15.21 . The dot-dashed diagonal line shows the one-to-one relation, 
which is followed nicely by confirmed RGB stars. It is reassuring that the M31 RGB distribution generally follows the RGB PDF and not 
the dwarf PDF. 



also impose a strict cut to eliminate a few very blue stars 
that are inconsistent with the color of M31's RGB (sec 
Gilbert et al. 2006). The remaining sample (530 stars) 
are shown as shaded histograms in Figure 0] We stress 
that although it may not be clear from any individual di- 
agnostic whether or not a star is a dwarf or a giant, the 
combination of the ten PDFs presents an unambiguous 
answer. 

As a point of interest, we note that some other studies 
of M31 have isolated Milky Way dwarfs from M31 RGB 
stars based solely on a velocity cut. Using our PDFs 
for the two minor axis fields Hll and aO, we calculate 
that such a cut will introduce a significant amount of 
Milky Way dwarf star contamination into the putative 
RGB sample. For example, invoking only a velocity cut 
of Wrad < —100 km s _1 to isolate M31 RGB stars from 
Milky Way dwarfs will lead to a contamination fraction of 
~15% in our data. Similarly, for a velocity cut of w ra( j < 
-150 km s" 1 cut, over 10% of all classified M31 RGB 
stars would actually be Milky Way dwarfs. For fields 
that sample the extended disk of M31 in the N-E or N- 



W quadrants, the contamination rates could be higher 
due to the rotation of M31's disk which would cause the 
velocity histogram to overlap with the Milky Way dwarf 
distribution. This would of course be offset to some de- 
gree (perhaps completely) depending on the distance of 
these fields from M31's center (i.e., the density of M31 
RGB stars is much higher along the major axis). 

In Table 2, we present the numbers of stars in each field 
that were determined to be secure M31 RGB stars and 
Milky Way dwarfs (measured using a (Li) > 0.5 cut) from 
the procedure outlined above. Further analysis from this 
point will involve only the cleaned M31 RGB sample. 

4. KINEMATICS OF SAMPLE 

In this section we briefly present kinematics of M31 
RGB stars in some of our fields and demonstrate that 
the data are consistent with a broad Gaussian centered 
at M31's systemic velocity. As can be seen from Fig- 
ure ^ nve °f our pointings are located on the minor 
axis (Hll, aO, m6, m8, and mil), two of the point- 
ings are located on the giant southern stream (H13s and 
a3), and five pointings are removed from the minor axis 
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TABLE 2 



Mask R No. Sci. No. M31 RGB No. Milky Way No. M31 RGB Comments 





(kpc) Targets 1 


o+ 2 
btars z 


Dwarfs^ 


Spheroid Stars 




Hll 


12 


277 


106 


18 


106 




rilos 


21 


272 


104 


20 


21 


giant southern stream field 


aO 


30 


268 


67 


30 


67 




a3 


33 


248 


68 


12 


20 


giant southern stream field 


al3 


60 


151 


18 


23 


18 




d3 


69 


242 


60 


89 


6 


And III dwarf spheroidal field 


al9 


81 


71 


4 


23 


4 




m6 


87 


147 


9 


19 


9 




bl5 


95 


139 


6 


25 


6 




m8 


121 


115 


1 


24 


1 




d2 


145 


280 


84 


38 





And II dwarf spheroidal field 


mil 


165 


295 


3 


82 


3 





1 Some targets were observed on multiple masks. 
2 Numbers given indicate secure M31 RGB stars ((Li) > 0.5) and 
Milky Way dwarfs ((Li) < -0.5) - see §|3]and Gilbert et al. (2006). 




Overall Likelihood: (L,) Overall Likelihood: <L,) 



Fig. 4. — The weighted average likelihood distinguishing M31 
RGB stars ((Li) > 0.5 — shaded) from Milky Way dwarfs for each 
of our twelve fields. The fields are arranged in order of increasing 
projected distance from M31. 



(al3, d3, al9, bl5, and d2). The data for the minor 
axis fields and the off axis pointings are likely fair rep- 
resentations of the bulge and halo of M31 (see 
However, the targeted stars in the giant southern stream 
(H13s and a3) and dwarf satellites (d2-And II and d3- 
And III) have very different kinem atics than hot bulg e 



or halo stars. As d e monstr ated i n llbata et al 



iGuhathakurta et "all l)2006bfl . and iKalirai et al 



(2004). 

2006), 



the stream consists of a kinematically cold population of 
stars blueshifted relative to M31 (in H13s stream stars 
have v ra d < —380 km s _1 , in a3 stream stars have w ra d < 
—400 km s" 1 ). The stream is well removed, both az- 
imuthally and radially, from all of our other pointings 
and does not contaminate these fields. Similarly, d2 and 



Fig. 5. — Radial velocity histograms for confirmed M31 RGB 
stars in six representative bulge and halo fields, arranged in order 
of increasing distance from M31. A scaled Gaussian with the char- 
acteristic velocity and dispersion of our most populous field, Hll, 
has been overlain on the data to show that the distribution of stars 
in each of these fields is roughly consistent with a broad Gaussian 
centered near M31's systemic velocity. 



d3 are dominated by stars belonging to the two dwarf 
satellites And II (-230 km s _1 < w rad < -150 km s _1 ) 
and And III (-390 km s" 1 < u rad < -325 km s" 1 ). For 
each of these four fields we easily remove the contribu- 
tion of the stream and satellite galaxies by making these 
velocity cuts, leaving a pure M31 RGB bulge and halo 
sample. 

Figure El presents the velocity histograms for six fields 
in M31 that are not contaminated by the giant south- 
ern stream or one of the dwarf satellites. Velocity his- 
to grams for the giant southern s tream can be found 
in IKalirai etTafl l)2006|) (H13s) and iGuhathakurta et alJ 
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(2006b) (a3). The fields shown here have been randomly 
chosen and sample the entire distance range over which 
we measure M31 RGB stars. The two minor axis fields 
with a large sample of stars, Hll and aO, clearly demon- 
strate that the majority of these stars form a broad, 
hot component centered near M31's systemic velocity. 
For our most populated field (Hll), we find u ra d ~ 
—330 km s _1 with a dispersion of a v ~ 90 km s _1 (skew 
= 1.3). The velocity histograms for our outermost fields 
are affected by small number statistics. However, al- 
most all confirmed M31 RGB stars in these fields (e.g., 
al9, m6, bl5, and mil) have kinematics consistent with 
M31's systemic velocity, w rac i ~ —300 km s _1 . To demon- 
strate this, we have overlaid a scaled Gaussian with the 
mean and dispersion of the Hll field in each of these 
panels (dashed curve). 

The combination of Figures 01 and [S] shows that we 
have detected a population of RGB stars belonging to 
M31's bulge and halo in eleven of our twelve fields (the 
exception being field d2, which is dominated by stars 
belonging to M31's dwarf satellite, And II). The final 
starcounts for these populations, in each field, are given 
in column 6 of Table 2 (this is the number after elimi- 
nating foreground dwarfs, M31 stream, and M31 satellite 
galaxy stars). Adding up all of the data, our sample con- 
sists of 261 bona fide M31 RGB bulge and halo stars, 47 
of which have R > 60 kpc. We stress that prior to this 
work, very few (if any) M31 halo RGB stars have been 
spectroscopically confirmed at these large radii. 

5. METALLICITY MEASUREMENTS 

We determine the metallicities of stars in our sample 
using two independent techniques, photometrically from 
the (I, V—I) CMD and spectroscopically from the Call 
triplet (A - 8500A). 

5.1. Photometric Metallicity Determination 

In Figure we present the CMD for our entire M31 
RGB bulge and halo sample. For this purpose, we 
first converted the Washington System (M, T%) pho- 
tometry into Johnson- Cousins (J, V— I ) ma gnitudes 
using the relations in IMaiewski et alJ (|2000|) . The 
systematic error in the slope of the color conversion 
from these relations i s less than 3.6%, as measured by 
IMaiewski et al.l l)2000|) . Also shown are several theoreti- 
cal isochrones ranging in metallicity from [Fe/H] = —2.31 
- +0.49, all with an age of 12 Gyr and [a/Fe] = 0.0 
ijVandenberg. Bergbusch. fe Dowlerl feOOS). These have 
been adjusted to a distance of 783 kpc ((m — M)q 
'2 1.17:. Most of our sample is indeed confined within the 
bounds of the most metal-rich and metal-poor isochrones 
on the CMD. For reference, the bolder isochrone has a 
metallicity of [Fe/H] = —1.0. The dashed line at the top 
of the isochrones indicates the tip of the RGB. 

We compute photometric metallicities ([Fe/H] p hot) for 
these stars by first measuring the length in the CMD 
of a segment that extends from the most metal-poor 
isochrone to the most metal-rich. This is done at 35 
points along the set of isochrones, extending from the 
base of our RGB sample in FigureElto well above the tip 
of the RGB (by linearly extrapolating the isochrones in 
brightness and color). The dashed curve on Figure [S] rep- 
resents one of these length segments. We then normalize 



20 - 




(v-D 



Fig. 6.— CMD of all M31 RGB bulge and halo stars. 
Most stars are confined within a grid o f theoretical isochrones 
IVandenberg. Bergbusch. fc D owlcr 20Q^) of age 12 Gyr spanning 
the metallicity range, [Fe/H] = —2.31 to +0.49 (bold isochrone: 
[Fe/H] = —1.0). A few outliers are found above the tip of the RGB 
(dashed line) . The arrows indicate the direction in which we ex- 
trapolated metallicity measurements for a few stars. Section 15.11 
discusses the determination of these metallicities for this sample. 



these 35 segments by their total length, hence comput- 
ing an index, X, equal to the fractional distance that a 
point is away from the most metal-poor isochrone. This 
X parameter is now a smooth function of metallicity, for 
a given Y range, the distance above the base of the RGB 
in our sample (extending from at the base to 1 at the 
tip of the RBG). Photometric metallicities are derived 
by measuring each star's X position and interpolating 
that within the relation that is appropriate given its Y 
value. Although we do this for many Y ranges (bins of 
size 0.1), we find that the X versus [Fe/H] relations are 
not very sensitive to this parameter and vary smoothly 
and slowly over the entire magnitude range of RGB stars 
in our sample. Metallicities for confirmed RGB stars that 
lie outside the range of the isochrones are derived by ex- 
trapolating the X versus [Fe/H] relations. As the CMD 
in Figure El shows, this mild extrapolation was required 
for only a few stars. 

We also tested our [Fe/H] p hot measurements derived 
above using two independent sets of isochrones. Inter- 
polating metallicitie s within a grid of the Padova models 
ijGirardi et al.l 12002]) and the Yale-Yonsei models (Y 2 - 
Demarque et al. 20 04) give consistent values to those 
derive d using the IVandenberg. Bergbusch. fc Dowlerl 
(2005) models. For a typical star, the mean [Fe/H] 
varies by less than 0.15 dex depending on which isochrone 
set is adopted. We find that the Padova models sys- 
tcmatical ly yield slightly more metal-rich values than 
both the IVjiid^nberg^F jcrgbusch. & Dowlcr (2005J) and 
IDemarau e et alJ (I2004D models. We choose to adopt 
the IVandenberg. Bergbusch. fc Dowlerl l)2005|) models as 
they also provide a grid of a-enhanced isochrones over a 
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Fig. 7. — Comparison of spectroscopic ([Fe/H] spcc ) versus pho- 
tometric ([Fe/H] p hot) metallicities of all M31 RGB bulge and halo 
stars. We have used a minimum bin size of 0.2 dex in [Fe/H] p hot to 
bin the sample, while ensuring a minimum of 20 stars in each bin. 
The data span a wide metallicity range and the two independently 
measured quantities are found to be in good agreement with one 
another over most of this range. For reference, a 1:1 relation is 
shown as a dashed line. 



Fig. 8.— The distribution of M31 RGB stars in the bulge 
(R < 20 kpc, open circles), crossover region (R ~ 30 kpc, 
crosses), and halo (R > 60 kpc, filled circles) on the CMD 
show differences. In 8 16.21 we provide evidence that these differ- 
ences are related to the MDF of stars in M31's different struc- 
tural components. Solid curves represen t theoretical isochrones 
IVandenberg. Bergbusch. fc Dowlerl 12003) with [Fe/H] = -2.31 
(left) and 0.0 (right). The dashed curves are discussed in 8 16.21 



broad metallicity range (see § 16 .2[) . 

5.2. Spectroscopic Metallicity Determination 

Independent of the method discussed above, we also 
determine metallicities for all of our M31 RGB bulge 
and halo stars using their spectra. This procedure re- 
lies on measuring the equivalent widths of the three 
Can absorption lines (see Figure EJ- The strengths 
of each of these lines are combined to yield a re- 
duced equ ivalent width according to the presc ription de- 
scribed in lRutledee. Hesser. fc Stetson! l)1997l) . This re- 
duced equivalent width is then calibrated empirically 
based on Galactic globular clust er RGB stars to yield 
[Fe/H] spec llRutledee et al.l 119971). Furthe r detai ls on 
this procedure are given i n iGilbert et alJ 1)200(11) and 
ICu hatha,kurta~et7hl BoOBbl) . 

In Figure [7| we present a comparison of the two inde- 
pendently determined metallicity measurements for our 
sample. Given the larger scatter in [Fe/H] spcc , we have 
binned our sample by using a minimum bin size of 0.2 
dex in [Fe/H] p hot while ensuring >20 stars in each bin. 
The individual [Fe/H] p hot and binned [Fe/H] spcc mea- 
surements are found to be in good agreement with one 
another over most of the metallicity range (1:1 line is 
shown), indicating that there are no systematic varia- 
tions in our metallicity scales. This ag rees with the re- 
sults of lE.eitzel fc G uhatha kurtal l)2002f > who also find a 
nice agreement between photometric and spectroscopic 
[Fe/H] measurements for RGB stars. 

We also note here that most classes of systematic 
biases/errors in our measurement of [Fe/H] p hot and 



[Fe/H] spec (e.g., age errors, systematic errors in the mod- 
els, etc.) will tend to affect different fields similarly. 
Therefore, although the absolute [Fe/H] p hot for any given 
star may have a large error (e.g., due to photometric 
calibration or incorrect distance modulus), the relative 
comparison of [Fe/H] p hot between different fields within 
our sample is more accurate. In S I7.1I we discuss several 
possible biases. 

6. ANALYSIS 

6.1. The Bulge and Halo Samples 

As Table 2 shows, our M31 coverage ranges from a field 
located at R = 12 kpc (Hll) to a field located at R — 
165 kpc (mil). Given this very large range in radius from 
M31's center, and our small number statistics in the out- 
ermost fields, we separate our fields into three broad re- 
gions. We define the bulge sample to be that represented 
by the Hll and H13s fields, the crossover region by the 
aO and a3 fields, and the halo by the al3, d3, al9, m6, 
bl5, m8, d2, and mil fields (d2 containing no halo mem- 
bers). This separation is justified naturally for several 
reasons. First, as discussed in §^ several authors have 
shown that there is no obse rved metallicity gradient in 
M31 out to R ~ 30 kpc (e.g.. lDurrell. Harris, fc PritchetJ 
2004). One study, Irwin et al. (2006), suggests that 
ther e is no metallicity gradient out to ^45 kpc . Sec- 
ond. iGuhathakurta et all <|2006a|) and llrwin etafl l|2005D 
have recently shown that a power-law R~ 2 profile domi- 
nates beyond this radius whereas a de Vaucouleurs i? 1 / 4 
surface brigh tness describes region s interior to this ra- 
dius. IGuhathakurta et al.l (|2006af) further demonstrate 
that the crossover of M31's bulge and halo occurs near 
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Fig. 9. — Normalized MDFs of confirmed M31 RGB bulge and halo stars in each of our twelve fields. A general trend indicating the 
presence of metal-poor stars in the halo of M31 is seen. As a guide, the dotted line is fixed at [Fe/H] = —1.0 in all panels. An arrow is 
indicated in the panel for m8 as we have re-n ormalized the MDF for this field. The cumulative distributions corresponding to these twelve 
fields are shown in Figure [Tfjl In Figure [TT1 we bin all of the halo data (R > 60 kpc) together to produce a more statistically significant 
sample from which we measure the mean metallicity and dispersion of M31's halo. Note: field d2 does not contain any M31 RGB stars 
that are not members of And II. 



this radius. Therefore, R ~ 30 kpc is a natural choice for 
separating the inner and outer samples. In addition to 
the two fields located near this intermediate radius (aO 
and a3), two of our best sampled fields are located inte- 
rior to this radius (Hll and H13s) while the remaining 
fields are located beyond this radius. 



6.1.1. Color- Magnitude Diagrams 

Figure M illustrates the CMD for M31 RGB stars in 
each of the three broad radial bins defined above. Also 
shown are two theoretical isochrones (solid curves - Van- 
denberg, Bergbusch, & Dowler 2005) with [Fe/H] = 
-2.31 (left) and [Fe/H] = 0.0 (right). As we showed in 
Figure El the distribution of M31 RGB stars nicely fol- 
lows the shapes of the isochrones. However, we can now 
see that some clear differences exist between stars be- 
longing to the inner (R < 20 kpc, open circles) , crossover 
(R ~ 30 kpc, crosses), and outer (R > 60 kpc, filled cir- 
cles) regions of M31. We show quantitatively below that 
the M31 RGBs in each of these regions are significantly 
different and that these differences relate to the MDF 
of each population. The dashed curves will be discussed 
later. 



6.2. Metallicity Distribution Function and Radial 
Trends 

Before presenting our MDFs for confirmed M31 RGB 
stars in each of the inner, crossover, and outer regions 
defined above, we construct MDFs for each of our twelve 
fields independently. These are shown in Figure [5] along 
with a fixed guide at [Fe/H] = —1.0 (dotted line). The in- 
nermost fields, Hll and H13s, contain mostly stars with 
intermediate/metal-rich compositions. The aO and a3 
fields at R ~ 30 kpc are dominated by a much broader 
MDF that extends to more metal-poor stars while still 
containing a metal-ric h component. This MDF n icely 
supports the finding in lGuhathakurta et all l)2006aj) that 
a roughly equal mix of bulge (metal rich) and halo (metal 
poor) M31 stars reside at R ~ 30 kpc (see §E2J>. Al- 
though individually affected by small number statistics, 
the outermost fields (second and third rows) in our study 
are dominated by a large number of metal-poor stars 
([Fe/H] phot < — 1)- Directly comparing these fields to the 
inner fields, Hll and H13s shows a very obvious metal- 
licity gradient in M31. Less than ^4% of all stars in the 
two inner fields (R < 21 kpc) combined have ([Fe/H] p hot 
< —1) whereas almost 65% of stars in the outer fields 
(R > 60 kpc) are this metal poor. Ostheimer (2003) 
found evidence for similar shifts in the proportions of 
metal rich and metal poor stars as a function of radius 
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Fig. 10. — Cumulative distributions of [Fe/H] p hot f° r M31 RGB bulge and halo stars in each of our twelve fields (solid curves). The 
dashed curves represent fixed Gaussians with [Fe/H] p hot = —1.3 and —0.4 to be used as guides. The same trend seen in the discrete 
distributions in Figure|5]is confirmed here. The outer fields are dominated by more metal-poor stars relative to the inner fields. Note: field 
d2 does not contain any M31 RGB stars that are not members of And II. 



from a purely photometric analysis of all giant star can- 
didates in these same fields. 

In Figure ITU1 we present cumulative [Fe/H] p hot distri- 
butions for each of the twelve fields. The distributions 
confirm the conclusions drawn above. The outermost 
halo fields contain metal-deficient stars relative to the 
inner fields. As a guide, we have also shown two Gaus- 
sian distributions with means at [Fe/H] p hot = — 1.3 and 
—0.4 and a — 0.4 (dashed curves). 

In Figure ^2 we present both the MDFs (top) and 
the cumulative Gaussians (bottom) for the bulge (Hll 
and H13s), crossover region (aO and a3), and the halo 
(al3, d3, al9, m6, bl5, m8, d2, and mil) of M31. All 
histograms have been normalized to unity for clarity. 
These data clearly show that we have indeed detected 
a metallicity gradient in M31. For R < 20 kpc, we find 
([Fe/H] phot ) = -0.47 ±0.03 (a = 0.39), for R ~ 30 kpc 
we find ([Fe/H] phot ) = -0.94 ± 0.06 (a = 0.60), and for 
R > 60 kpc, we find ([Fe/H] phot ) = -1.26 ± 0.10 (a = 
0.72). Gaussian models with the mean and dispersion of 
the stars comprising the MDF in each of these groups are 
overlaid on the differential histograms as dashed curves. 
For completeness, we note that our spectroscopic metal- 
licity for the outer halo (([Fe/H] spcc ) = -1.24 ± 0.12, 
a = 0.85) is in excellent agreement with the measured 
photometric metallicity. 

The CMDs for stars belonging to each of these three 
regions were presented in Figure El To illustrate the dif- 



ferences in the three CMDs, we have overlain isochrones 
in each panel of Figure |S1 with the approximate mean 
metallicity of the respective population as determined 
above (dashed curves). 

The photometric metallicities given above have been 
derived assuming isochrones with [a/Fe] = 0.0. If we 
assume that M31's halo is a-enhanced, with [a/Fe] = 
+0.3, the mean metallicity of the outer halo is found to 
be ([Fe/H] phot ) = -1.48 ±0.11 (a = 0.73). 

6.3. Kolmogorov-Smirnov Tests 

To test whether or not the MDFs of the inner, 
crossover, and outer regions of M31 (as defined above) 
truly differ significantly from one another, we have 
applied the two-s ided Kolmogorov-Smirnov (K-S) test 
l)Press et al.111992^1 . The K-S test makes no assumption 
about the distribution of the data and is therefore insen- 
sitive to possible spurious biases from arbitrary binning 
of data. The test looks at the cumulative fraction of 
each histogram and returns a statistic, D, which mea- 
sures the maximum vertical deviation between the two 
curves. This statistic yields a significance level prob- 
ability P, that the null hypothesis (i.e., that the two 
data sets are drawn from the same distribution) is true. 
Therefore, small values of P show that the two MDFs are 
significantly different whereas larger values (e.g., >0.1) 
indicate that the two MDFs may have been drawn from 
the same distribution. 
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Fig. 11. — (Top) — Normalized MDFs for each of the bulge (Hll and H13s), crossover region (aO and a3), and halo (al3, d3, al9, m6, 
bl5, m8, d2, and mil) regions. The dashed curves represent Gaussians with parameters fixed to the measured mean and dispersion of 
each MDF. The vertical dotted line is held fixed at [Fe/H] = —1.0 to guide the eye. (Bottom) — Cumulative distributions, corresponding 
to the above panels. The dashed curves represent fixed Gaussians with mean [Fe/H] = —1.3 and —0.4 (a = 0.4). The data unequivocally 
show a metallicity gradient in M31. 



We find that the R < 20 kpc inner M31 MDF is sig- 
nificantly different from the R ~ 30 kpc MDF, yielding 
a tiny probability P — 9.9xl0 -8 of the two being drawn 
from the same parent distribution. We also find that 
the R ~ 30 kpc MDF is significantly different from the 
halo MDF, R > 60 kpc. For this, the two-sided K-S 
test returns a probability of P = 7.5xl0 -3 . Additional 
tests confirm that the aO and a3 MDFs that have been 
combined to form the R ~ 30 kpc population are likely 
drawn from the same distribution, P = 0.23 (at least 
it can not be shown that they are drawn from different 
distributions). Directly comparing the inner fields, Hll 
and H13s, with any of the outer fields al3, d3, al9, m6, 
bl5, m8, or mil shows them to be significantly different. 
Therefore, the K-S tests confirm our earlier conclusion 
that we have detected a metallicity gradient in M31. 

In Figure ^] we present the radial metallicity distri- 
bution for M31. The mean photometric metallicity in 
each of our fields is plotted as a function of radius (open 
circles). The mean metallicity of the combined data for 
the inner, crossover, and outer regions is also plotted as 
larger filled circles. The data clearly show a trend of de- 
creasing metallicity as the distance from M31 increases. 
For comparison, we also show the mean spectroscopic 
metallicity in each of our inner, crossover, and outer sam- 
ples as filled squares. These measurements confirm that 
M31's outer halo is dominated by stars much more metal- 



poor than the bulge. Table 3 summarizes our final MDF 
results. 

7. DISCUSSION 

7.1. Systematic Errors and Measurement Biases 

As discussed in §[IJ the detection of a radial metallicity 
gradient in M31 can directly constrain the global nature 
of chemical enrichment during the galaxy formation pro- 
cess. Our finding, however, is in contrast to previous 
observational work on M31 which has found no signif- 
icant radial [Fe/H] gradient within R < 3 kpc (e.g.. 
Dun-e ll et al. 2001, 2004). Even the recent llrwin et alJ 
( 2005) photometric survey found no significant color gra- 
dient for M31 RGB candidates out to R ~ 55 kpc, but 
uncertainties in statistical background subtraction make 
it difficult to interpret the significance of this null result. 
Given the relevance of a metallicity gradient and the ap- 
parent discrepancy with earlier studies, we now address 
and quantify some potential biases in our study. Under- 
standing these effects is crucial to determining whether 
or not a spurious metallicity gradient could have been 
produced by our target selection and/or analysis proce- 
dure. 

The target selection for the two innermost fields, Hll 
and H13s, is different from the outer fields in the sense 
that DD051 filter pre-selection was not used. Although 
any selection effect from this may alter the overall shape 
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TABLE 3 



Field 


{R) 1 


No. M31 RGB 


{[Fe/H] 


phot) cr 




(kpc) 


Spheroid Stars 






Bulge 


14 (12-21) 


127 


-0.47±0.03 ( 


-0.57±0.04) 0.39 (0.43) 


Crossover 


31 (30-33) 


87 


-0.94±0.06 ( 


-1.12±0.07) 0.60 (0.68) 


Halo 


81 (60-165) 


47 


-1.26±0.10 ( 


-1.48±0.11) 0.72 (0.73) 



1 Weighted average of the projected radial distance based on the 
numbers of stars in each individual field that were grouped to- 
gether, and the range of those individual distances (in parenthe- 
ses). 

2 Photomctric mctallicities and dispersions calculated assuming 
[o/Fe] = 0.0. Values in parenthesis assume [a/Fe] = +0.3 
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Fig. 12. — The mean metallicity of M31 RGB stars in each of 
our fields as a function of increasing radius (open circles). The 
error bar on the most metal-poor data point, from the m8 field, 
has been omitted as it is based on only one stars and therefore 
is not reliable. The mean metallicity in the three coarser bins 
representing the inner, crossover, and outer regions of M31 are 
shown as larger filled circles. Spectroscopic metallicity measure- 
ments for these three components are shown as filled squares (the 
most metal-poor point strongly overlaps the photometric metallic- 
ity point). The data clearly shows that the halo of M31 is more 
metal-poor than the inner regions. 



of the [Fe/H] distribution, there should not be any rel- 
ative bias among different DD05Tselected fields or, for 
that matter, among non-I?Z?051-selected fields. To in- 
vestigate further whether a global bias exists between 
the inner and outer fields we can take advantage of 
the fact that two of our fields target the giant south- 
ern stream of M31. The target selection for one of 
these fields, H13s, does not involve DD051 photome- 
try, whereas for the other field, a 3 , it d oes. It is reas- 
suring t hat both the iKalirai et all 112 006) study of H13s 
and the IGuhathak urta et alJ l|2006bf) study of a3 found 
a substantial population of high metallicity stars in the 
giant southern stream despite the I < 22.5 magnitude 
cut used for spectroscopic target selection in both cases 



(see details below) . A detailed comparison between these 
two fields can reveal whether DD051 preselection biases 
the sample against the most metal-rich stars. 

In Figure ^5] we present the MDFs in H13s and a3, 
including only bona fide stream members. These have 
been selected by imposing velocity cuts on the data. In 
H13s, Kali rai et al.l l|2006|) show that the stream consists 
of stars moving wit h v r!tt \ < —460 km s _1 . Similarly, 
IGuhathakurta et alJ (J2006b) show that at the position 
of a3 stream stars have w ra d < —425 km s -1 . There 
is some evidence for secondary velocity peaks in both 
of these data sets that may or may not be related to 
the stream, but we ignore these components here (see 
Kalirai et al. 2006 for further information). The result- 
ing two MDFs are remarkably similar and show a strong 
peak of metal-rich stars with a small tail to metal-poor 
stars. For H13s, this MDF is well represented by a double 
Gaussian (shown as a dashed curve in the upper panel of 
Figure lT3*|l . We now superimpose this same model on the 
a3 data (dashed curve in lower panel), and find that the 
model that fits H13s needs to be adjusted only slightly 
by A[Fe/H] p hot = +0.1 dex to achieve the best fit of the 
stream MDF in a3: this was determined by fixing the rel- 
ative spacing, standard deviations, and relative fraction 
of stars in the two components of the double Gaussian 
and only fitting for an overall offset. This shifted double 
Gaussian is shown as a dotted curve in the lower panel. 
The comparison of these two (differently selected) stream 
samples shows that any systematic [Fe/H] offset between 
them is very small. Since there are several other sources 
of systematic errors in [Fe/H] at this level or greater, we 
ignore this 0.1 dex offset; in fact, this apparent [Fe/H] 
offset may well be a result of the fact that our a3 sample 
probes deeper down the RGB luminosity function than 
the H13s sample (see discussion below). We note that 
if we did apply this potential offset between DD051- 
selected and non-DD051-selected data, the radial metal- 
licity gradient in M31 would become stronger. 

We note that an intrinsic metallicity gradient in the 
giant southern stream is unlikely over such a small ra- 
dial extent (R — 21-33 kpc). Several investigations have 
shown that the stream is young and was produced less 
than an orbital period ago (e.g., < 2 Gyr — Ibata et al. 
2004). The radial extent above represents a very small 
difference in dynamical age and it would be difficult for a 
metallicity gradient to have developed within this time- 
frame. 

Several additional biases exist in our data that would 
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tend to weaken an observed gradient compared to any 
true radially-outward decrease in the mean metallicity of 
M31: 

• The tip of the RGB is fainter for more metal- 
rich populations (see Figure EJ). Thus, our mag- 
nitude cut for spectroscopic target selection (20 < 
Iq < 22.5) can bias the sample against the inclu- 
sion of the most metal-rich RGB stars (Reitzel & 
Guhathakurta 2002). This bias, if present, is ex- 
pected to be stronger for the bulge than for the 
halo. This is because our target selection procedure 
(see Guhathakurta et al. 2006b) includes some- 
what fainter targets, on average, in the outer fields 
due to the sparseness of luminous M31 RGB stars 
in these regions. 

• The same PDF (see §|2Jl is being used to separate 
M31 RGB stars from foreground Milky Way dwarfs 
in all fields. This selection would tend to pick out 
RGB stars from the same region of the [Fe/H] p hot 
(or location in the CMD or V— I color) space at all 
radii (see Gilbert et al. 2006 for details). Given 
that our RGB PDF is based on "training set" stars 
drawn mostly from the inner fields, these PDFs 
favor higher [Fe/H] p hot values. In fact, Figure|3](d) 
clearly shows that the RGB training set PDF does 
not sample metal-poor stars very well. Therefore, 
the true radial trends in metallicity may be even 
larger than observed. 

• We have assumed a constant age of 12 Gyr for al l 
stars in M31. However, the iBrown et at] l|2003l) 
study indicates that the bulge contains a substan- 
tial fraction of intermediate age stars. By contrast, 
one expects the halo population to be old by anal- 
ogy with the Milky Way's halo. If in fact the inner 
region of M31 probed in this study is younger on 
average than the halo, then our photometric metal- 
licities have been underestimated (i.e., assigned to 
be more metal-poor than they truly are) in this re- 
gion. In fact, the comparison between [Fe/H] p hot 
and [Fe/H] sp cc in Figure ITSl suggests this. There- 
fore the adoption of younger isochrones in the bulge 
would lead to higher metallicities for these stars 
and a larger metallicity gradient. A shift in age 
from 12 to 6 Gyr translates to a +0.3 dex offset in 
[Fe/H] phot . 

• We have assumed that all stars in M31 are not 
enhanced in their a-element abundances, ([a/Fe]). 
Stars in the Milky Way's halo are known to be a- 
enhanced ([a/Fe] = +0.3) and it is generally be- 
lieved that this is a result of the stars in the halos 
of galaxies forming in early "bursts". If the M31 
halo is also a-enhanced, then our use of [a/Fe] 
= isochrones has overestimated the metallicity 
for stars in this component (i.e., assigned to be 
more metal- rich than they truly are). We deter- 
mined earlier that the mean metallicity of our M31 
RGB halo sample is in fact ~0.2 dex more metal- 
poor if we assume [a/Fc] = +0.3 (i.e., we calculate 
([Fe/H] pho t) = -1.5 for M31's halo under this as- 
sumption). 
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Fig. 13. — MDFs for giant southern stream stars in H13s (top) 
and a3 (bottom). The two MDFs are found to be very similar. 
We find that the best fit double-Gaussian model to H13s (dashed 
curve) needs to be shifted only by A[Fe/H] phot = +0.1 dex to fit 
the distribution of stars in the a3 field. This offset may result from 
a bias related to the way obj ects were selected in the inner and 
outer fields in our study (see S I7.1l for more information). 



These additional biases strengthen our detection of a 
metallicity gradient in M31 are suggest that, if anything, 
the detected gradient is likely smaller than the true gra- 
dient. 

7.2. Bulge/Halo Ratio 

The analysis presented in this paper suggests that, over 
a large range of radial distances (12-165 kpc), there is 
a metallicity gradient within M31. This brings up the 
question: Is this metallicity gradient intrinsic to the 
bulge and/or halo, or are we seeing two overlapping M31 
structural components that are each homogenous but dis- 
tinct from each other in terms of their chemical abun- 
dance properties? One way to test these two scenarios 
is to consider the crossover region at 30 kpc, where we 
find a mean metallic i ty of ( [Fe/H] p h t,) = —0.94 ± 0.06. 
iGuhathakurta et alJ (|2006afl showed that at this radius 
there is roughly an equal mix of bulge and halo stars 
based entirely on fits to the minor-axis surface bright- 
ness profile of M31. If our data are consistent with the 
latter scenario (e.g., a homogenous bulge and halo), then 
the 30 kpc MDF should in fact be an equal superposi- 
tion of the distinct bulge and halo MDFs. To test this, 
we randomly selected 47 bulge stars and grouped them 
together with the same number of stars in the halo sam- 
ple. The resulting composite MDF is compared to the 
30 kpc MDF derived from the aO and a3 fields and found 
to be very similar. A two-sided K-S test (see S I6.3|I be- 
tween this composite "bulge+halo" MDF and the 30 kpc 
MDF returns a probability of the null hypothesis being 
true of P > 0.50 (i.e., of the two samples being drawn 
from the same parent population). 

We therefore conclude that our data are consistent 
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with a chemically homogenous bulge and halo whose rela- 
tive contributions change systematically with radius and 
give rise to the observed metallicity gradient. This is not 
to say that M31's bulge and halo are perfectly homoge- 
nous in a large-scale radial sense — only that our present 
samples lack the size and precision to detect any subtle 
metallicity gradients that may be intrinsic to the bulge 
or halo. For example, Figure IT51 shows a very mild trend 
that an inverted metallicity gradient may be present in 
the halo of M31. However, splitting our 60-165 kpc 
halo fields into two bins (e.g., R < 100 kpc and R < 
100 kpc) suggests that such a trend is statistically not 
significant (marginally) . This measurement is difficult to 
make given the small number of stars. Future samples 
will no doubt be able to rectify this situation and allow 
us to characterize radial trends in the intrinsic chemical 
abundance properties of the bulge and halo in detail. 

7.3. Relation of our Fields to M31 's Disk 

None of the fields used in this study are contami- 
nated by M31's extende d disk to any significant de- 
gree. Ilbata et alJ l|2005f) find that this disk is a low- 
surface brightess, kinematically cold (velocity dispersion, 
^30 km s" 1 ) structure at R = 15-40 kpc and is an ex- 
tention of the inner disk. 

Our fields tend to lie close to M31's minor axis, and are 
therefore well removed from the major axis. The most 
susceptible fields in our study that could potentially suf- 
fer from this contamination are the two innermost point- 
ings, Hll and H13s. The projected distance of H13s is 
21 kpc, which in the plane of the disk is >80 kpc. There- 
fore tfils_iieki_2^ocated well beyond the extendend disk 
and lKalirai et alJ l)2006j) have already ruled out any sig- 
nificant contribution to this field by a smooth disk pop- 
ulation. We now focus on Hll, which is located at R — 
12 kpc in the bulge of M31 (i.e., 50 kpc in the disk). We 
can rule out disk star contamination in this field based 
on several independent arguments, 

• IGuhathakurta et "all l)2006aj) present the surface 
brightness profile of M31 and show that at the po- 
sition of Hll, disk stars are expected to be greatly 
outnumbered by bulge stars. The bulge contri- 
bution to Hll is at least an order of magnitude 
greateer than the disk contribution. 

• Figure [S] (Rich et al. 2006, in preparation) shows 
that the velocity histogram of Hll is well repre- 
sented by a broad Gaussian centered near M31's 
systemic velocity. There is no evidence for a 
30 km s _1 kinematically cold population (as seen 
in the fields that probe the extended disk of M31 in 
Ibata et al. 2005). In fact, the velocity dispersion 
of M31 RGB stars in Hll is measured to be ap- 
proximately three times higher than that of M31's 
extended disk. A maximum likelihood analysis of 
the Hll radial velocity data indicates that it is im- 
possible to hide a substantial fraction of the stars in 
a kinematically cold (30 km s _1 ) component (Re- 
itzel et al. 2006, in preparation). Although this 
extended disk population dominates Ibata et al.'s 
RGB samples in fields near the major axis, their 
two minor axis fields, F05 and F07, also show no 
evidence of this population, as noted by the au- 
thors. 



• Another scaling argument based on kinematics can 
be used to rule out the presence of a substantial 
fraction of M31 disk member s in the Hll field. Fig- 
ure 9 in iKalirai et al.l l)2006|) shows radial velocity 
histograms for Hll (R = 12 kpc, S-E minor-axis), 
H13s (R = 21 kpc, giant southern stream), and 
H13d (R = 25 k p c, N- E major-axis). Consistent 
with Ilbata et all l)2005|) . our H13d field shows a 
kinematically cold population of stars, presumably 
M31's disk. This field lies at 25 kpc in the disk 
plane, and also at 25 kpc in the bulge (as it is on the 
major axis). We find that two-thirds of the RGB 
stars in this field belong to the cold component 
(disk) and one-third are part of a dynamically-hot 
(bulge) component (Reitzel et al. 2006, in prepa- 
ration). By comparison, our Hll field is located at 
50 kpc in the disk plane and 20 kpc in the bulge (ac- 
counting for the 5:3 flattening of the bulge). Thus, 
the Hll field is twice as far out in the disk as H13d 
(ten versus five disk scale-lengths) and at a slightly 
smaller radial distance in the bulge. Moving out in 
the disk by five exponential scale lengths reduces 
the disk contribution by over two orders of magni- 
tude. Thus, very few, if any, disk members should 
be present in our Hll field. 

• T. Brown et al. (2006, in prep) have recently 
found evidence for younger stars in M31's disk 
that are not seen in Hll. By obtaining ultra-deep 
HST / ACS observations of stars in both of these 
fields (directly overlapping our pointings), Brown 
et al. are able to reach the main-sequence turnoff 
in M31's bulge and disk. The reconstructed star 
formation histories of the these two components in- 
dicate some clear differences. For example, the disk 
field contains a younger main sequence and is gen- 
erally more metal-rich than the bulge field. There- 
fore, these independent observations also suggest 
that our Hll field is not contaminated with disk 
stars and rather represents M31's bulge. 

Therefore, the surface brightness, kinematics, ages, 
and metallicities of stars in Hll argue against them be- 
ing part of a disk like component in M31. We note that 
our next closest field, aO, is located at Rsim 30 kpc in 
M31's halo (i.e., over 130 kpc in the disk) and therefore 
d oes not sample the ex tended disk. 

iWorthev et alJ l|200 51 make the radical hypothesis that 
a thick disk dominates the entire region within R < 
50 kpc in M31. This hypothesis was put forward to ex- 
plain the high mean metallicity seen in previous studies 
of M31, but was not based on any kinematical data. Re- 
itzel et al. (2006, in preparation) show that the stellar 
kinematics in the inner regions of M31 are inconsistent 
with Worthey et al.'s hypothesis, and confirm that the 
Hll and H13s samples used in this work represent M31's 
bulge population instead. 

7.4. Metallicity Distribution of the M31 Bulge: 
Comparisons to Other Studies 

Although the halo sample of M31 RGB stars pre- 
sented in this work is unique, a limited number of stud- 
ies have probed the MDF of the bulge of M31 (although 
these studies have often referred to this component as 
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the "halo"). The best s tudy comes from ultra-deep 
HST/ ACS observations bv lBrown etaD £2006) who have 
studied both the bulge and giant southern stream. These 
data reach faint enough to detect the lower RGB, hor- 
izontal branch, and main-sequence turnoff of M31. At 
these fainter magnitudes, there is minimal foreground 
dwarf star contamination. The HST CMDs of these 
two fields show striking similarities over all of the above 
phases of stellar evolution suggesting that the ages and 
metallicities of stars in the two fields are very similar. 
Furthermore, our Hll and H13s Keck/DEIMOS spectro- 
scopic pointings have been carefully chose n to directly 
overla p these two HST /ACS fields. In H rown et all 
( 2003|, it was shown that the Hll CMD contains both 
metal-poor and metal-rich stars ([Fe/H] > —0.5). Given 
the striking similarity between the HST CMDs in these 
two fields, the H13s field must have a comparable MDF. 
Our measurement of the MDF in Hll and H13s confirm 
this (see Figure 1X0)1 . both showing a distribution of stars 
near [Fe/H] 0.5. 

The MDF of M31's bulge has also been probe d pho- 
tometrically by iDurrell. Harris, fc Pritcheti (I2001D (R ~ 
20 kpc) and IDurrell. Harris. & Pritchetl l)2004|) {R ~ 
30 kpc). These measurements are based on wide field 
CFHT data. Durrell et al. compute the MDF by inter- 
polating the location of stars on the CMD within a grid 
of a-enhanced RGB models. This is done both for the 
target fields and for well-removed control fields. The final 
MDF is obtained by subtracting the control field MDF 
from the target field MDF. The results, in both fields, in- 
dicate that M31's bulge is dominated by relatively high 

metallicity stars, ([m/H]) 0.5 (([Fe/H]} = -0.8, for 

[a/Fe] = +0.3). As Table 3 shows, our inner field is 
located interior to Durrell et al.'s 20 kpc pointing, and 
our crossover field is located slightly exterior to Durrell 
et al.'s 30 kpc pointing. Our measured metallicities, as- 
suming the same [a/Fe] = +0.3 adopted by Durrell et 
al., are found to nicely bracket Durrell et al's metallic- 
ity (see Table 3). This suggests that Durrell et al. have 
succeeded in statistically eliminating Milky Way dwarf 
stars with their use of an appropriate control field. Of 
course, this becomes more difficult as the distance from 
M31 increase s , and therefore the methods discussed in 
IGilbert et all <|2006j) become more important. 

Other studies of M31 (e.g., Worthey et al. 2005, 
Bellazzini et al. 2003) have similarly concluded that 
the bulge is dominated by metal-rich stars (([Fe/H]) ~ 
—0.5). Both these studies, and those discussed above, 
vary from one another in terms of both the source of the 
data and the types of models used to determine metal- 
licity. Therefore, external comparisons of these results 
to our MDFs for Hll and H13s may not yield a perfect 
agreement and are not as powerful as relative compar- 
isons of the MDF at various radii from a single data set 
(as presented in this paper). However, it is encourag- 
ing that such an external comparison does in fact yield 
similar metallicities and thus suggests that net effect of 
systematic errors in our MDFs is small. 

7.5. The M31 Halo Versus the Milky Way Halo 

In §H we discussed the assembly of massive galaxies 
such as the Milky Way and M31 in the context of hi- 
erarchical merging of smaller galaxies. Several predic- 
tions arise from this theoretical framework that have, 



until now, been challenged by our understanding of M31. 
Both structurally and in terms of chemical enrichment, 
the M31 "halo" was believed to be different from our 
own Galactic halo and different from the predictions of 
generally-accepted halo formation model s. The wide- 
field imaging survey of lOstheimerl l|2003D was the first 
to find preliminary evidence of a metal-poor stellar halo 
in M31: it used DD051 selection to greatly reduce fore- 
ground contamination by Milky Way dwarf stars and was 
thereby able to probe out to larger radii than previous 
studies. By coupling the Ostheimer DDQ51 photometry 
with K eck/DEIMOS spectroscopy, iGuhatha kurta et "all 
( 2006aJ) showed that M31 does possess a stellar compo- 
nent that is structurally distinct from its bulge and whose 
radial extent (R > 160 kpc) and surface brightness pro- 
file (~ R~ 2 ) resemble those of our Galaxy's stellar halo. 

In this paper we have taken another step toward 
bridging the apparent disparity between M31 and the 
Galaxy/galaxy formation models by showing that this 
newly-discovered stellar halo in M31 is chemically dis- 
tinct from its bulge, and is in fact quite metal poor: 
([Fe/Hlnhnt.) = -1-5 ± 0.1 (<t = 0.7). For comparison, 
iMorrison et alJ l)2003|) have recently measured metallici- 
ties from spectra of Milky Way halo RGB stars located 
at distances between 15 and 83 kpc from the Galactic 
center. The mean metallicity of this sample is measured 
to be ([Fe/H]) = —1.6 (<r = 0.6), in good agreement with 
earlier studies. Therefore we find that the mean metal- 
licity, and spread, of the M31 halo is similar to that of 
the Galactic halo. 

The results presented in this paper show that a survey 
out to 20-30 kpc is not sufficient to properly isolate the 
M31 halo. In the Milky Way, the halo begins to dominate 
the thin disk at a radius interior to 10 kpc, locally (Ma- 
jewski 1993). Therefore, our sampling distance in M31 is 
not comparable to the Milky Way. M31 likely contains a 
substantially larger bulge than the Milky Way and this 
bulge dominates over the halo in all previous studies of 
M31 that were limited to the central few degrees of the 
galaxy. Our data show clear evidence of a halo and bulge 
in M31 that are distinct from each other in terms of both 
structure and chemical abundance properties. 

Prior to the detection of 47 M31 halo stars in this anal- 
ysis, very few (if any) spectroscopically confirmed RGB 
stars had been detected in M31 with R > 60 kpc. Fu- 
ture observations will allow us to gather larger bulge and 
halo samples to establish whether or not these compo- 
nents contain intrinsic radial metallicity gradients. A 
detailed study of the true halo of M31 will require new 
observational strategies to gain significant insights into 
this important, but tenuous/elusive, stellar population. 

8. CONCLUSIONS 

Imaging and multiobject spectroscopic surveys of M31 
have recently brought to light several unanticipated re- 
sults. By isolating M31 RGB stars from Milky Way fore- 
ground dwarf contamination using a new technique, we 
are able to probe further into the outskirts of M31 than 
previous studies (R — 12-165 kpc). In a separate pa- 
per we established the existence of a true stellar halo 
in M31, thereby resolving a long-standing concern that 
the surface brightness profile of M31 looks different from 
the Milky Way. In this paper, we analyze the metallic- 
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ity distribution function of stars in twelve fields span- 
ning this large radial extent. Our data show a radial 
metallicity gradient in the spheroid (bulge and halo) of 
M31. The outer region of M31 (R > 60 kpc), domi- 
nated by its halo, is composed of stars deficient in met- 
als relative to those in the inner region. Based on the 
sample of 47 stars, we find a mean chemical abundance 
of ([Fc/H] phot = -1.48 ± 0.11 (a = 0.73) for the newly 
discovered halo of M31 (assuming [ct/Fc] = +0.3). This 
confirms predictions from current galaxy formation mod- 
els which suggest that the inner regions of large galaxies 
should be chemically enriched relative to the outer re- 
gions. 
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